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摘 要 
甘于 经 典 动力 学 模型 和 模型 参数 自动 优化 算法 , 本 文 分 析 了 全 国 当前 累计 确诊 感染 病例 数 大 于 100 
的 24 个 省 市 自治 区 , 以 及 湖北 省 除 神农 架 以 外 的 16 个 地 市 从 2020 年 1 


(截止 到 2020 年 2 月 16 日) 

月 20 日 至 2 月 16 日 的 累计 确诊 病例 数 , 并 对 相 应 省 市 地 区 疫情 可 能 的 结束 时 间 和 总 感染 人 数 的 进行 
了 长 期 预测 . 我 们 的 研究 表明 ， 在 目前 严防 严 控 措施 下 , 全 国 大 部 分 省 市 的 疫情 将 于 2 月 底 前 基本 结束 ， 
而 湖北 省 内 疫情 也 有 望 于 3 月 中 旬 结 束 , 但 是 武汉 市 的 疫情 可 能 要 持续 到 4 月 初 . 通过 对 比 公开 数据 
与 预测 值 , 我 们 建议 加 强 对 黑龙 江 、 河 北 、 江 西 、 安 徽 、 贵 州 和 四 川 六 省 , 以 及 湖北 省 内 武汉 、 荆 州 、 
鄂州 、 随 州 、 天 门 和 恩施 六 个 地 市 的 的 监控 , L 防疫 情 死 灰 复 燃 . 此 外 ， 分 析 结果 提示 ， 在 疫情 发 展 前 
期 ， 天 津 、 河 北 、 重 庆 、 四 川 、 海南 和 广西 等 省 , 以 及 湖北 省 下 辖 多 个 地 市 可 能 存在 聚集 性 感染 ， 这 
人 于 疫 后 进一步 的 流行 病 学 调查 确认 . 关键 字 : 新 型 冠状 病毒 肺炎 ”种群 动力 学 ”动力 学 模型 ”疫情 
抽 测 

1 引言 


2019 年 底 , 新 型 冠状 病毒 肺炎 (COVID-19, 以 下 简称 新 冠 肺炎 ) 疫情 在 武汉 先 爆发 ， 其 后 迅速 昌 延 到 全 
司 力 至 海外 地 区 . 截止 到 论文 完成 之 日 ， 我 们 仍然 身 处 严重 的 疫情 之 中 .除了 传 纪 充 的 医学 和 生物 学 手段 ， 
基于 统计 或 数学 建 模 的 方法 能 够 帮助 我 们 有 效 地 分 析 并 预测 疫情 的 基本 流行 特征 、 拐 点 、 结 束 时 间 以 及 
感染 人 数 , 从 而 为 控制 疫情 传播 , 制定 相关 的 防 控 政策 法 规 提 供 理论 支持 . 
疫情 早期 的 流行 病 学 研究 主要 集中 于 对 基本 再 生 数 [1], 传染 速率 [2], 倍增 时 间 [3] 等 关键 性 参数 的 获 
取 和 估计 . 根据 这 些 参数 , 并 结合 传染 病 演化 动力 学 模型 ， 如 SI、 SIR、 SEIR 等 , 就 可 以 对 疫情 基本 传播 


情况 进行 定量 的 分 析 和 预测 [4, 5]. 本 文中 , 我 们 报告 了 基于 推广 的 经 典 SEIR 模型 [6], 对 全 国 当 前 累计 

确诊 感染 病例 数 大 于 100 (截止 到 2020 年 2 月 16 日 的 24 个 省 市 自治 区 , 以 及 湖北 省 除 神农 架 以 外 

的 16 个 地 市 从 2020 年 1 月 20 日 至 2 月 16 日 的 累计 确诊 病例 数 (来 自 国家 及 相关 省 市 卫生 健康 委员 
会 官 万 公 数据 ) 的 分 析 结 果 , 以 及 对 相应 省 市 地 区 疫情 可 能 的 结束 时 间 和 总 感染 人 数 的 长 期 预测 . 

我 们 所 使 用 的 差分 方程 模型 如 图 1 所 示 , 其 中 包含 了 易 感 者 (5S)、 不 易 感 者 (有 效 防护 , P)、 潜 伏 

者 (处 于 潜伏 期 暂 不 具备 传染 能 力 , 媚 、 感 染 者 (具备 传染 能 力 且 尚未 被 隔离 , 1) 和 确诊 者 (确诊 且 被 隔 

离 , @) 五 类 基本 人 群 . 模型 中 的 四 个 转移 速率 参数 {a, 6,y-1,6-!} 均 为 非 负 的 , 分 别 表示 防护 率 、 感 染 


率 、 平 均 潜伏 期 (根据 [7], 范围 为 0~24 天 , 中 位 数 为 3 天 ; 另 参考 [3], 疫情 前 期 95% 置信 区 间 为 4.1~7.0 
天 ) 和 平均 隔离 响应 时 间 (根据 [8,9], 全 国平 均 为 4.95 天 ). 其 中 参数 a 的 引入 , 主要 是 考虑 到 国民 健康 意 
识 的 提高 (如 佩戴 口音 ) 和 防 控 措施 的 加 强 (追踪 密切 接触 者 和 封闭 社区 ) 等 ， 使 得 易 感 人 群 逐 渐 减 少 , 这 
也是 与 经 典 SEIR 模 型 的 主要 差别 之 一 . 文献 [6] 中 所 考虑 的 治愈 和 死亡 人 数 , 为 了 简单 起 见 , 本 文中 不 作 


根据 文献 [7] 报道 , 模型 中 平均 潜伏 期 固定 为 中 位 数 3 和 0 -1 = 3), 除 潜伏 者 和 感染 者 人 数 外 ， 
模型 中 其 余 变 量 初 值 均 可 从 公开 数据 中 直接 获得 . 这 样 模 型 er {QB,6-!} 和 待定 隐 变 量 初 始 


* 通 讯 作者 邮件 地 址 : zcamhl@tsinghua.edu.cn( 洪 柳 ), penglr.whQ@qq.com( 萌 梁 榕 ) 


易 感 者 (9 | 个 >| 潜伏 者 (BD) | a 加 于 者 有 各 
a Si+1) =50)— SOW — si), 
不 易 感 者 (P) ] | a D ] Eli+1)= ECG Pe —yE(D), 
| I(i+1)=10)+yEG)— 610), 


确诊 者 (0) od+D= (0) + 610). 


Figure 1: COVID-19 传 染 模 型 示意 ( 左 ) 及 对 应 差分 方程 ( 右 ). 


值 (0), 7T(0)， 将 通过 极 小 化 模型 预测 值 与 已 知 时 间 序 列 Q(t) 之 间 的 累计 相对 误差 来 获得 并 优化 , 而 优 
化 问题 的 全 局 解 将 通过 模拟 退火 法 目 动 求解 . 参数 和 初始 值 对 模型 预测 的 影响 可 参见 [6]. 


2 ”全国 主要 省 市 疫情 基本 情况 分 析 和 预测 


我 们 首先 研究 了 全 国 当 前 累计 确诊 病例 数 大 于 100 (截止 到 2020 年 2 月 16 日 ) 的 24 个 省 市 自治 
基本 疫情 情况 . 分 析 表 明 , 全 国 各 省 市 普遍 的 隔离 响应 时 间 在 4~5 天 , 这 与 文献 [8,9] 中 报道 的 国内 患 
从 发 病 到 确诊 的 平均 时 间 是 4.95 天 相符 . 我 们 注意 到 东北 三 省 、 西 南 四 省 市 ( 除 西 藏 ) 和 西北 地 区 陕 
的 隔离 时 间 普 遍 较 长 ( 见 附录 参数 表 ), 特别 是 黑龙 江 省 、 云 南 省 、 四 川 省 和 陕西 省 , 这 可 能 对 疫情 防 


村 
+ 
cE 了 3 [Xl 
信号 孙 


来 一 定 挑战 . 拟 合 得 到 的 感染 率 参数 普遍 较 高 , 反映 了 新 冠 肺炎 极 强 的 感染 性 (如 日 本 的 钻石 公主 号 邮轮 ， 
在 未 充分 隔离 密闭 空间 中 的 感染 率 目 前 达到 了 621/3711 16.7%%, 参见 [10]). 具有 显著 统计 性 差异 的 省 
市 包括 天 津 、 河 北 、 重 庆 、 四 川 、 海南 和 广西 , 其 中 天 津 市 的 确诊 案例 主要 集中 在 宝 城区 百货 大 楼 附近 
x 域 (截止 2 月 20 日 占 比 为 55/130 s 42.3%, 见 [11]). 据 此 我 们 怀疑 上 述 省 市 , 特别 是 海南 省 和 重庆 市 
可 能 也 存在 类 似 情况 , 还 需要 进一步 通过 流行 病 学 调查 确认 . 

此 外 预测 结果 表明 , 湖北 之 外 的 大 部 分 省 市 的 疫情 拐点 (基本 再 生 数 小 于 1) 集中 于 1 月 29 日 到 2 
月 1 日 ( 见 图 2), 稍 早 于 基于 新 增 确诊 病例 的 2 月 4 日 全 国 拐点 [12]. 其 有 省 的 拐点 定 于 1 月 26 日 ， 
EL 于 该 省 内 较 少 的 、 散 发 性 病例 数 , 而 河北 和 贵州 的 拐点 在 2 月 六 后 , 一 定 程度 上 反映 了 这 两 


人 
3 日 以 后 ， 

份 早期 疫情 的 严峻 性 . 模型 预测 大 部 分 省 市 ee 小 于 5 和 小 于 1 确 
定 ) 集 中 于 2 月 中 下 旬 ( 见 图 1), 而 河北 、 河 南 、 江 7 `、 广东、 重庆 、 云 南 、 贵州 的 疫情 结束 
日 期 则 需 推迟 1 周 左右 (一 般 在 3 月 上 旬 ), 特别 是 四 川 和 黑龙 江 i 省 的 疫情 持续 时 间 会 更 长 , 需要 进 
步 加 强 监 控 . 上 述 结果 与 我 们 之 前 对 全 国 的 总 体 预 测 是 一 致 的 [6]. 最 后 , 对 比 预测 的 总 感染 和 人数 和 当前 报 
告 数据 (并 选择 预测 的 可 能 感染 人 数 高 于 目前 报道 人 数 10% 的 省 市 ), 我 们 建议 加 强 对 黑龙 江 、 河 北 、 江 

、 安 徽 和 贵州 5 个 省 的 疫情 监控 . 全 国 ( 除 湖北 省 外 ) 最 终 可 能 感染 人 数 在 1.2~1.5 万 . 


3 ”湖北 省 下 辖 地 市 疫情 基本 情况 分 析 和 预测 


于 目前 疫情 主要 集中 于 武汉 市 和 湖北 省 内 地 区 , 我 们 进一步 分 析 了 湖北 省 内 除 神农 架 以 外 的 16 个 地 市 
的 基本 疫情 情况 , 并 对 未 来 疫情 可 能 走向 做 了 初步 预测 ( 见 图 3). 首先 由 于 2 月 12 日 - 14 日 临床 确诊 病例 
数据 的 加 入 , 使 得 湖北 省 内 大 部 分 地 区 累计 确诊 病例 有 了 大 幅度 增加 , 这 也 对 模型 拟 合 和 参数 估计 提出 了 
更 高 的 挑战 . 综合 考虑 各 种 因素 , 我 们 较为 可 靠 地 确定 了 相关 参数 和 初始 值 , 并 发 现 相对 省 外 , 虽然 大 部 分 
湖北 省 内 地 市 都 有 大 量 的 感染 病例 告 ， 同时 却 有 更 多 地 市 的 感染 率 参 数 表现 出 显著 的 统计 性 差异 , 其 

包括 鄂州 、 研 门 、 咸 宁 、 和 十堰 、 天 门 、 恩 施 和 洪江 ( 见 附录 参数 表 ). 这 意味 着 上 述 地 区 之 前 可 能 存在 着 聚 
集 性 感染 . 其 次 ， 从 隔离 响 应 来 看 , 武汉 、 荆 州 、 鄂 州 、 随 州 、 天 门 和 恩施 等 地 的 3 均 隔离 响应 时 间 都 超 
过 6 天 , 远 远 高 于 全 国平 均 水 平 ， 而 前 面 四 个 城市 恰恰 是 前 一 阶段 疫情 传播 最 为 迅速 的 地 区 之 一 . 天 门 和 
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9a oa eR ya 6 和 人 Sg sy amd 
大 陆地 区 华北 华南 西北 西南 
城市 北京 | 天 津 | 河北 | 山西 | 广东 | 广西 | 海南 | 陕西 | 四 川 | 重庆 | 云南 | 贵州 
预测 拐点 日 期 1.29 5 2.4 1.30 1.30 1.29 2.1 1.29 1.29 1.28 1.26 2.4 
预测 结束 日 期 2.22-2.29 | 2.25-3.4 | 2.28-3.6 | 2.20-2.28 | 2.28-3.5 | 2.20-2.27 | 2.25-3.5 | 2.24-3.6 | 3.5-3.17 | 2.28-3.8 | 2.24-3.8 | 2.29-3.10 
预测 最 终 感染 人 数 | 370-450 | 130-160 | 320-390 | 130-160 | 1300-1600 | 220-270 | 180-220 | 230-280 | 490-600 | 540-650 | 160-200 | 170-200 
大 陆地 区 华东 东北 华中 
城市 山东 | 安徽 | 江苏 | 上 海 | 浙江 | 福建 | 黑龙江 | 辽宁 | 吉林 | 河南 | 江西 | 湖南 
预测 拐点 日 期 1.31 1.31 2 1.28 1.28 1.29 2 1.28 2.1 2.1 2 1.30 
预测 结束 日 期 226340 32300 | 313.80 22332 0 227350 332333 | 3.10-3.21 | 2.18-2.27 | 2.23-3.3 | 2.293.6° 323.10 | 22532 
预测 最 终 感染 人 数 | 530-650 | 990-1200 | 630-770 | 320-390 | 1100-1400 | 270-330 | 530-650 | 120-140 | 99-120 | 1200-1500 | 950-1200 | 960-1200 
未 预测 内 蒙古 、 甘 肃 、 青 海 、 新 疆 、 宁 夏 、 西 藏 3 单列 | 湖北 
. AN 确 六 二 -Fx 、 
Figure 2: 全 国 24 省 市 公布 累计 确诊 人 数 、 预 测 确诊 人 数 及 预测 潜在 感染 者 人 数 的 演化 趋势. ! 所 列 为 
和 /从 二 本 
分 省 市 的 疫情 拐点 、 结 束 期 (上 限 / 小 于 5 人 ， 下 限 小 3 -1 人 ) 和 总 感染 人 数 预 测 数据 总 结 . 
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城市 武汉 孝感 黄冈 莉 州 鄂州 随州 惠阳 黄石 
预测 拐点 日 期 2.9 al 1.30 2.1 2.3 1.30 1.29 1.31 
预测 结束 日 期 4.10-4.19 3.9-3.17 | 3.3-3.10 | 3.17-3.28 | 3.17-3.28 | 3.3-3.12 | 2.27-3.5 | 2.29-3.7 
预测 最 终 感染 人 数 ”| 46000-57000 | 3300-4000 | 2700-3300 | 1600-2000 | 1300-1600 | 1200-1500 | 1100-1300 | 960-1200 
城市 宜昌 荆门 咸宁 十 堰 仙桃 天 门 恩施 潜 江 
预测 拐点 日 期 1.30 1.31 2.3 1.31 2.2 2.4 1.25 2.3 
预测 结束 日 期 3.1-3.10 3.8-3.16 | 3.8-3.16 | 2.29-3.7 3.1-3.9 3.19-4.1 3.9-3.24 | 3.5-3.14 
预测 最 终 感 染 人 数 870-1100 | 950-1200 | 860-1100 | 610-740 | 550-670 | 500-620 | 250-310 180-220 
未 预测 神农 架 


Figure 3: 湖北 省 下 辖 16 个 地 市 ( 除 神农 架 ) 公 布 累 计 确诊 人 数 、 预 测 确诊 人 数 及 预测 潜在 感染 者 人 数 的 演 
化 趋势 . 表 中 所 列 为 各 地 市 的 疫情 拐点 、 结 束 日 期 (上 限 小 于 5 人 ， 下 限 小 于 1 人 ) 和 总 感染 人 数 预测 数据 总 
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隔离 响应 时 间 可 能 源 于 当地 较为 偏僻 的 地 理 位 置 , 建议 进一步 加 强 防 控 响应 . 此 儿 
4 隔离 响 9 时 间 ,经 接近 全 国平 均 水 平 , 这 表明 在 上 述 地 区 由 于 前 段 时 间 积 
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Forecasting the ending of COVID-19 epidemics by generalized 
SEIR models 


Changjing Zhuge & Dongyan Zhang & Wuyue Yang & Liangrong Peng & Liu Hong 


Abstract Here we report the analysis of epidemic data from Jan. 20th to Feb. 16th, 2020 in 24 provinces 
in China，whose total infected cases are larger than 100 till 02/16/2020, as well as 16 cities in Hubei province 
(the most severely affected area) except Shennongjia, based on dynamical models and automatical algorithms for 
parameter optimization. We forecast the COVID-19 epidemics in most provinces in China will end up soon before 
February 29th, while those for Hubei province (except Wuhan city) will be closed by the middle of March. The 
epidemic in Wuhan will continue to the beginning of April. And we suggest further close attentions should be paid 
to six provides, including Heilongjiang, Hebei, Jiangxi, Anhui, Guizhou and Sichuan, as well as six different cities, 
including Wuhan, Jingzhou, Ezhou, Suizhou, Tianmen and Enshi, in Hubei province. Moreover, it is hinted that 
clustering infection might be happened in Tianjin, Hebei, Chongqing, Sichuan, Hainan and Guangxi provides, 
and many cities inside Hubei province during the spreading of COVID-19, which needs further validation by 
epidemiological investigations in the future. 


Keywords COVID-19, population dynamics, dynamical equations, epidemic forecast 


城市 常 驻 人 口 总 数 ( 根 


城市 据 2018 年 湖北 省 人 口 防护 率 | 感染 率 潜伏 局 离世 初始 暴 初始 感 预测 拐点 预测 结束 日 期 项 测 最 终 

年 鉴 ， 单 位 ， 万 人 ) C pb yl Or 露 人 Bo | 染 人 万 日 期 感染 人 数 
武汉 1108 0.1 1 3 6.25 1423 637 Feb. 9 Apr.10-Apr.19 51513 
孝感 492 0.17 1 3 5.4 388 138 Feb. 1 Mar. 9-Mar. 17 3632 
黄 闷 633 0.18 0.97 3 4.7 313 200 Jan. 30 Mar.3-Mar.10 2974 
荆州 559 0.18 0.92 3 7.6 158 106 Feb. 1 Mar.17-Mar.28 1800 
鄂州 108 0.14 0.63 3 7.2 189 75 Feb. 3 Mar.17-Mar.28 1459 
随州 222 0.26 1 3 6 262 100 Jan. 30 Mar.3-Mar.12 1355 
训 阳 557 0.28 0.89 3 5 176 185 Jan. 29 Feb. 27-Mar.5 1210 
黄石 247 0.21 1 3 5 104 92 Jan. 31 Feb. 29-Mar.7 1070 
宜昌 414 0.21 1 3 5.7 106 75 Jan. 30 Mar.1-Mar.10 965 
莉 门 290 0.1 0.41 3 5.8 208 74 Jan. 31 Mar.8-Mar.16 1053 
咸宁 254 0.1 0.52 3 5.3 76 90 Feb. 3 Mar.8-Mar.16 960 
十 堰 341 0.12 0.59 3 4.6 79 58 Jan. 31 Feb. 29-Mar.7 673 
仙桃 114 0.15 1 3 5.2 47 25 Feb. 2 Mar.1-Mar.9 611 
天 门 127 0.1 0.42 3 8.6 72 38 Feb. 4 Mar.19-Apr.1 560 
恩施 338 0.12 0.12 3 10 100 67 Jan. 25 Mar.9-Mar.24 279 
潜 江 97 0.1 0.56 3 6 14 22 Feb. 3 Mar.5-Mar.14 204 
神农 架 未 预测 


大 陆地 区 | 城市 0 Sk 人 Be < 0 约 *% i a "0 预测 结束 日 期 0 
东北 黑龙 江 37.7 0.158 1 3 7.3 49 17 Feb. 3 Mar. 10-Mar. 21 589 
辽宁 43.6 0.21 0.93 3 6 21 8 Jan. 28 Feb. 18-Feb. 27 128 
吉林 27 0.138 1 3 5.7 7 3 Feb. 1 Feb. 23-Mar. 3 110 
北京 21.5 0.158 1 3 4.5 26 26 Jan. 29 Feb. 22-Feb. 29 406 
华北 天 津 15.6 0.086 0.73 3 4.3 8 6 Feb. 2 Feb. 25-Mar. 4 148 
河北 75.6 0.084 | 0.745 3 3.9 24 13 Feb. 4 Feb. 28-Mar. 6 358 
山西 37.2 0.158 1 3 5.1 10 8 Jan. 30 Feb. 20-Feb.28 145 
内 蒙古 未 预测 
机 河南 96.1 0.151 1 3 4.4 118 63 Feb. 1 Feb. 29-Mar. 6 1346 
江西 46.5 0.134 1 3 4.9 60 41 Feb. 1 Mar. 2-Mar. 10 1058 
湖南 69 0.149 1 3 3.8 115 48 Jan. 30 Feb. 25-Mar. 2 1063 
湖北 单列 
山东 100.5 0.13 1 3 4.2 48 19 Jan. 31 Feb. 26-Mar. 4 593 
华东 安徽 63.2 0.146 1 3 4.8 62 57 Jan. 31 Mar 2-Mar. 9 1100 
江苏 80.5 0.132 1 3 4.9 36 27 Feb. 1 Mar. 1-Mar. 8 701 
上 海 24.2 0.202 1 3 5.5 37 25 Jan. 28 Feb. 23-Mar. 2 351 
浙江 57.4 0.206 1 3 5.1 190 72 Jan. 28 Feb. 27-Mar. 5 1255 
福建 39.4 0.225 1 3 5.0 26 33 Jan. 29 Feb. 23-Mar. 3 299 
广东 113.5 0.134 1 3 4 130 49 Jan. 30 Feb. 28-Mar. 5 1424 
华南 E 49.3 0.13 0.848 3 4 30 11 Jan. 29 Feb. 20-Feb. 27 249.3 
海南 9.3 0.089 0.541 3 5.1 12 17 Feb. 1 Feb. 25-Mar. 5 196 
Ll 83.4 0.201 | 0.817 3 7.7 68 44 Jan. 29 Mar. 5-Mar. 17 543 
西南 重庆 31 0.144 | 0.595 3 5.8 107 39 Jan. 28 Feb. 28-Mar. 8 595 
云南 48 0.316 1 3 8.2 37 20 Jan. 26 Feb. 24-Mar. 8 182 
贵州 36 0.12 1 3 5.9 9 3 Feb. 4 Feb. 29-Mar. 10 186 
西藏 未 预测 
陕西 38.6 0.266 1 3 7.2 38 25 Jan. 29 Feb. 24-Mar. 6 257 
西北 甘肃 
a 
本 未 预测 


